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I. INTRODUCTION 



The external gravitational field of a black hole, formed by collapse, settles down to the 
Kerr-Newman solution of General Relativity. These generic stationary black hole solutions 
have a rather curious feature, the inner horizon acts as the boundary of predictability for 
external initial value problems. However, Penrose |T| pointed out that this inner (Cauchy) 
horizon is a surface of infinite blueshift, and thus partial absorption of the radiative tail of 
the gravitational collapse leads to a divergent flux of energy along it. This suggests that the 
backreaction of perturbations on the geometry may generate unbounded curvature along the 
Cauchy horizon (CH). Detailed calculations, which analysed the evolution of gravitational 
perturbations confirmed the divergence of the energy flux |^, and for charged spherical 
black holes Poisson and Israel(PI) showed that a scalar curvature singularity does form 
at the CH when the above influx is accompanied by the outflux emitted from the collapsing 
star. Since they found the mass function diverges, they called it a mass inflation singularity. 
An exact solution, based on the work of PI was used by Ori to examine the nature of the 
singularity in some detail and he showed that an observer, who falls into the black hole, 
experiences finite tidal distortion at this singularity. 

The presence of a cosmological constant (of arbitrary magnitude) modifies this picture 
in a non-trivial manner, and serves to emphasize an aspect of the singularity along the 
CH which has so far been ignored. The leading divergence of the curvature scalars is not 
contained in the Weyl scalar, but rather in Raji^sR'^^^^ ■ 

For a generic perturbation of a Reissner-Nordstrom black hole in de Sitter space, the CH 
can be stable for some values of the charge (e) and the mass (m) [|],^. These configurations 
are quite rare though, and require \e\ > m. In p it was speculated that this spacetime 
might provide a counter-example to the conjecture that non-scalar curvature singularities 
are catalyzed into stronger scalar curvature singularities, once perturbed. This was based 
on the observation that the Weyl curvature of the solution to the Einstein field equations 
(for simple spherical perturbations) is dominated by the proper time integral of the energy 
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flux along the horizon. It was shown in P] that this flux behaves like T{v) ~ e^'^'^' and 



where w ^ oo on the CH. The constants and are the surface gravities of the inner 
and cosmological horizons respectively. Thus one finds values of the charge and the external 
mass of the black hole for which < K[ < 2kc, and the observed energy density grows 
without bound, but \&2 is finite. 

Although these observations are correct, a scalar curvature singularity does form along 
the CH contrary to the speculation in 0. We will demonstrate this fact using a model in 
which perturbations of a Reissner-Nordstrom-de Sitter black hole are modelled by cross- 
flowing streams of null dust, as in f^. The singularity is characterized by the divergence of 
RaisR'^^, where the dominant behavior is proportional to the influx of blueshifted radiation. 
We wish to emphasize that this behavior should be present in all spherically symmetric 
solutions, where the outflux from the star is continuous, and does not require a non-zero 
cosmological constant. However, the presence of a cosmological term shows clearly that 
the divergence of the mass function is not necessary for the presence of a scalar curvature 
singularity. 



where dfl"^ is the line element on the unit two-sphere, and latin indices a,b, . . . range over 
(0,1). Along with the scalar function r{x"') we introduce m{x°') , f [x"") and /t(a;") defined by 




(1.1) 



II. AN APPROXIMATE SOLUTION 



The spherically symmetric line element can be written as 



ds^ = da"^ + r'^dVt^, da'^ = gabdx°'dx^ , 



(2.1) 



g^'r^ar,, = f{x-) = 1 



2m(a;") A 



(2.2) 




K = dr f- 
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(2.3) 
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It is also convenient to use null coordinates U, V on the "radial" two spaces so that da'^ = 
—2e'^^dUdV and v = ^{U, V). It is now easy to generalize the derivation of the field equations 
in to include a cosmo logical constant. The result is 

r-ab + i^Qab = -4:7rrTab , (2.4) 

1 3e2 A 



□z/ = -d^K = — ^ I 2m - — + -r^ , (2.5) 
j-i y J- 3 / 

□m = -(47r)Vr,,T'^\ (2.6) 



where a semi-colon indicates the covariant derivative associated with da"^, and Tnh is the 2x2 



sub-matrix of the stress-energy tensor. 

For crossflowing null dust, the stress-energy tensor is 

where = —d^V and = —d^U are radial null vectors pointing inwards and outwards 
respectively, and, pj„ and pout represent the energy densities of the inward and outward 
fluxes. Covariant conservation requires that they are given by 

Let us imagine that the inflow is turned on at a flnite advanced time V = Vq and 
the outflow is turned on at f/ = Uq. In the pure inflow(outflow) regime the solution is 
an ingoing(outgoing) Vaidya-Reissner-Nordstrom-de Sitter spacetime with mass function 
m{V){m{U)). The global structure of the spacetime with cross- flow is shown in the Penrose 
diagram of Fig. 1. The loci of the three apparent horizons, in the pure inflow region, are 
given by the positive roots of / = 0. The ingoing sheet of the inner horizon is r = ri, and of 
the cosmological horizon is r = Tc. The outer apparent horizon settles down to the constant 
radius r = Ve, for asymptotically large v, and the function k evaluated at a particular horizon 
gives the surface gravity of that horizon. 

We choose V such that V = at the Cauchy horizon. Note that this does not specify 
the coordinate completely, but ensures that it may be regularly related to the Kruskalized 



advanced time associated with the inner apparent horizon. Our eventual aim is to construct 
an approximate solution which is valid as \^ ^ 0. 

Taking the trace of ( p.4| ) we find that r{U, V) also satisfies 

□r = -2K. (2.9) 

Using (|2.2| ), ( p.5|) and (|2^ ) we obtain the wave equation 

ailnire^")) = ^(Se' - - Ar^) , (2.10) 
which can be integrated formally to 

re^- = ag,{U)g2{V) exp [F{U, V)] , (2.11) 

with 

The functions gi{U) and (72(^) are determined by the boundary conditions along the null 
rays V = Vq and U = Uq, and a is a constant with the dimension of length. In |^ and [Q] it 
was shown that null coordinates exist such that g2{V) is well behaved as — 0, in particular 
we can rescale the coordinates so that (72(0) = 1. This guarantees that V = —(const.) x e~'^'''", 
where ki is the surface gravity of the static portion of the inner horizon. Similarly, provided 
we limit ourselves to values of U near to Uq we may set gi{U) = 1 without loss of generality. 

In order to proceed, we must estimate the behavior of the integral (|2.12|) . The crucial 
observation (which was made in but not carried to its full consequences) is that the 
integrand is negative for r < n, and we expect the leading behavior to come from near 
the CH. Thus if it diverges, the integral must diverge to negative values and re^" 0. As 
an immediate consequence of this, we see from ( |2.12| ) that this would require that r go to 
zero. Since we expect r{U,V) to be a well behaved function, with the slow contraction of 
ingoing null rays governed primarily by the outfiux from the star, we conclude that a good 
approximation to the metric coefficient is 



a exp {F{U, 0)] 



(2.13) 



near to = 0. 

With this ansatz we can rewrite 



as 



(2.14) 
(2.15) 



where F is treated as a function of U alone from here on. The solution of these equations 
is then given by 



/ ^ 



r 

IVo 









V-2 f^e' 




JUo 



dU' . (2.16) 



We still have not determined in detail what the function F{U, 0) is, but it is unimportant 
in the ensuing analysis. 

Before proceeding to analyse the curvature singularity which forms due to the infinitely 
blueshifted influx, we present the solution of equation (|2.6| ) for the mass function. In these 
coordinates and using (|2.13|) it is easy to integrate this to 



''^e-^'Lout(t/')l dU' r[LUV')]dV' 

JVo 



Uo 



+m(f/o, Vo)-(3 f [V'L,^,{V')] dV + 7 / [U'L,^,{U')] dU' , (2.17) 

JVn JUo 



''\v'L,^,{v')]dv' +7 n 

'Vb JUo 

where (3 and 7 are constants. The contribution from the pure inflow (outflow) region is finite 
for realistic perturbations, and it is only when both fluxes are present that the first term 
can lead to a divergence of the mass function. 

The functional form of the luminosities Lout{U) and Li^iV) is as yet unspecified, so 
we now specialize to the case of interest. In ^ and it was shown that for a generic 
perturbation the initial conditions on the influx should be 



L{V) ~ K{V)[-V/a] 



(2.18) 



near to the CH. This is proportional to the energy density of the influx, as measured by 
a free falling observer crossing the CH in the pure inflow region. Notice that it diverges 
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provided > Kc- The behavior of the prefactor in (p.l8|) is unimportant, so we take 



K{V) — > constant a.s V ^ (if the cosmological constant is zero, then = and K{V) 
is given by the analysis of Price as K{V) ~ | ln(— V"/a;)|~^). Since it is slowly varying 
compared to the other factor in ( |2.18| ), we can write 

J dVL{V) ~ K{V)[-V/a]^^^^. (2.19) 

The mass function is proportional to ( p.l9| ), so it will diverge as unless 2kc > ki. 



III. THE CAUCHY HORIZON SINGULARITY 



Our solution, given by (|2.13| ), ( |2.16|) and (|2.17|) , gives conditions under which the mass 



function is finite, and shows that the metric is well behaved provided r is bounded away from 
zero. In this section we examine the nature of the singularity that is present for a general 
outflux. There are three essentially different possibilities; (i) complete stability when Kc > i^u 
(ii) divergent influx but finite mass function when 2kc > > and (iii) both divergent 
influx and mass function when < n\. Defining a,=e^/m\ and h = Amf/3, where mi is 
the asymptotic value of the mass function in the pure inflow region, we have plotted these 
conditions in Fig. 2. Some useful relations involving the surface gravities are also given in 
appendix A. 

We wish to focus our analysis on case (ii) where no mass inflation has occurred despite 
the presence of the infinite energy density. At first sight, one is tempted to conclude that the 
spacetime is regular, since we know that the mass function characterizes the Weyl curvature 
scalar \E'2 oc m/r^ within spherical symmetry. However with a moment's thought one realizes 
that Eqn. ( p.7|) implies that RapR""^ diverges on the CH. Thus it should come as no surprise 
that the Kretchsmann invariant, given by 

\r^),uu-Fu{r^),u 



+ ... (3.1) 
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to leading order, also diverges like the measured energy density provided that the outflux 
from the star is non-vanishing. This shows that a scalar curvature singularity forms along 
the CH whenever a blueshifted influx and continuous outflux from the star are present. 



IV. CONCLUSION 

For completeness, we will take the Ori model limit of our solution. This is achieved by 
replacing the continuous outflux by a delta function source at U = Uq and then integrating 
over U where possible. The result is 
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rl-2fdV' I dV"L,^t{V") + f e'^'[V-2b]dU' , (4.1) 
m{V) ~ const, x / Vut(^')] dV - (3 / V'^out(V^')] dV + . . . (4.2) 
where h = , and is in complete agreement with M. Since Lout{U) = for f/ > f/g, it 

U=Uo 



is clear that the leading term in (|3.2| ) is absent, and the divergence of the curvature scalars 
is proportional to the mass function. Therefore this model misses the leading behavior 
discussed above, but we wish to stress that the tidal forces (and hence tidal distortion) 
experienced by a free-falling observer in both models are almost identical! 

In our opinion the result ( ^.2| ) is rather important. It shows that Lont{U) must be 
continuous in order to see the full nature of the singularity which forms along the CH. 
Based on this observation it is tempting to speculate that asphericities may act to change 
the nature of the singularity by pushing divergences (proportional to Lin(V")) into the Weyl 
curvature. Indeed, asymptotic and perturbative analyses of the singularity inside a more 
realistic black hole, suggest that the leading terms are in the radiative part of the Weyl 
curvature. It seems to us very important, then, to further investigate the internal structure 
of black holes in aspherical models. 

At the present time, the indications are that the singularity which forms inside a realistic 
black hole is null. A question which deserves some attention is, how generic is this null 
picture? Within spherical symmetry, we see that the singularity becomes spacelike when 

8 



r{U,V) 0. The decrease in r{U,V) (near the CH) is governed primarily by the outflux 
from the star [Eqn. ( |2.16| )], and thus the null portion of the singularity can be quite large. 
It is possible (likely) that the presence of shear in more realistic models will change this, 
provoking a spacelike singularity which is asymptotically null near P in Fig. 1. This would be 



consistent with the perturbative analysis in and, as pointed out by Yurtsever ||TD[ , would 
be similar to the situation in plane wave spacetimes where perturbations do not capture 
the generic structure for colliding waves. Finally we suggest that it is rather premature to 
discuss the strength of the singularity which forms, or the possible fate of an astronaut who 
falls into a black hole, since it is likely that tidal forces will be substantially enhanced if the 
singularity in more realistic models is spacelike. 

Some of these issues are currently under investigation and details will be presented else- 
where. 



ACKNOWLEDGMENTS 

It is a great pleasure to thank Werner Israel for fruitful comments and encouragement, 
and also Eric Poisson for constructive criticism which helped us clarify the exposition. D. 
N. is grateful for the awards given by External Affairs and International Trade Canada, 
Government of Canada Awards, administrated by the International Council for Canadian 
Studies. He also thanks DGAPA, UNAM for support. This work was also partly supported 
by the Natural Sciences and Engineering Research Council of Canada. Figures 2, (i) and 
(ii), and some calculations in this work were done using Mathematica© . 



APPENDIX A: THE SURFACE GRAVITIES 

In order to obtain the physical regions in the ab-plane in figure 2, and to examine the 
conditions on the surface gravities we found the following useful. 

Lake [Til] introduced a rescaled quantity x = r/mi, where mi is the mass of the black 



hole, so that the scalar function / is given by 
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m^x"^ bob 



The surface gravity at any horizon r is similarly given by 

1 



(Al) 



rriibx^ 



(A2) 



By construction the roots of (|A1|) satisfy < Xj < < Xc, and we know that the surface 
gravities are always greater than zero. These facts and ([A2|) tell us that 



Xi < 



3- 



< Xc < 



3 + yg^^ 



(A3) 



We now observe that for a = 9/8, we have Xe = 3/2 and since this satisfies ( |A1| ), b = 2/27. 



In fact for these values of a and b we have Xi = Xg = Xc |TT 



Finally, the coalescence of the two or more roots occurs when Kr = 0. Using ( |A1| ) gives 
, y'^-2y + a 3 ± V9 - 8a 



1/4 ' ^ 2 ' 

which allows us to plot the boundary of the regions of interest. 



(A4) 
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FIGURES 

FIG. 1. The global structure of the spacetime, with crossflowing null dust. EH is the event 

horizon, CH is the Caucliy liorizon and the cosmological horizon is at r = Tc- The lines U = Uq 
and F = Vb at which the fluxes are switched on are also shown. 

FIG. 2. The ab-plane. Along curve 1 Vc = r^, and along curve 2 ri = re- (i)The region of 
physical interest is bounded by the axes and the curves 1 and 2. (ii) A close up of the physical 
region. In the shaded portion no mass inflation occurs, but energy densities are divergent on the 
CH. Above this, both the mass function and the energy density diverge, and below is the region 
of complete stability 
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